Aging and reduced bulk conductance in thin films of the topological insulator Bi2Se3 
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We report on the effect of exposure to atmospheric conditions on the THz conductivity of thin 
films of the topological insulator Bi2Se3. We find: 1) two contributions of mobile charge carriers 
to the THz conductivity immediately after growth, and 2) the spectral weight of the smaller of 
these decays significantly over a period of several days as the film is exposed to ambient conditions, 
while the other remains relatively constant. We associate the former with a bulk response, and 
the latter with the surface. The surface response exhibits the expected robustness of the carriers 
from 2D topological surface states. We find no evidence for a third spectral feature derived from 
topologically trivial surface states. 



One of the outstanding issues in the field of topolog- 
ical insulators^^ has been the unique identification of 
the topological metallic surface states in transport mea- 
surements made on both single crystals and thin films. 
Signatures of the surface response have been contami- 
nated by appreciable bulk conduction. Only recently 
have clear signatures of 2D behavior been identified in 
transport experiments^^. I n particular, 2D behavior has 
been demonstrated using thin films of Bi2Se3 of differ- 
ent thicknesses and studying the thickness dependence 
of the transport properties®. This study showed that a 
contribution to the conductance was independent of the 
film's thickness, and thus 2D in origin. The bulk carriers 
dominated the conductance response for film thicknesses 
larger than 300 nm. These same films have been used to 
study the Kerr rotation induced by magnetic field from 
the surface states^. Even though the films were exposed 
for several days to atmospheric conditions, the Kerr re- 
sponse changed little, indicating the robustness of these 
states against aging due to exposure to atmospheric con- 
ditions. 

Recently, experiments that explored the effects of ex- 
perimental conditions on the response of the surface and 
bullf^ have been reported. Bianchi et al.^ studied the 
co-existence of trivial (i.e. without topological proper- 
ties) and topological surface states in single crystals of 
Bi2Se3 using photoemission experiments. During these 
experiments the sample was kept in ultrahigh vacuum, 
and it was argued that the trivial states are formed due to 
quantum confinement of the 3D electronic structure near 
the surface, due to so-called band-bending. In addition, 
King et alM^ showed that the largest Fermi momentum 
crossing of the surface states, and accordingly the total 
carrier density, changed with time on a timescale of a few 
hours, even as the sample was kept in ultrahigh vacuum. 
They argued that this change was due to adsorption of 
residual gases present in the ultrahigh vacuum chamber, 
such as CO, on the surface. Similar results were obtained 
by deliberate dosing the surface of freshly cleaved single 
crystals with water molecules by Benia et al.^. In this 
work, it was explicitly demonstrated that upon further 
water exposure, the Fermi energy increased up to values 



of « 0.6 eV, while at the same time topologically trivial 
surface states were induced due to confinement. These 
studies raise the question of what the effect is on the 
transport properties of topological insulators of exposing 
the them to atmospheric gases. 

In this letter we explore the effects of exposure to at- 
mospheric conditions on the surface states in thin films 
of the topological insulator Bi2Se3. We use time domain 
terahertz spectroscopy (TDTS) to study the finite fre- 
quency properties of thin films of Bi 2 Sc 3 grown by molec- 
ular beam epitaxy on sapphire (AI2O3) (0001) substrates. 
By using a spectroscopic probe, we are able to separate 
the effects of aging on the scattering rate (mobility) and 
on the spectral weight (the ratio of carrier concentration 
to the effective mass). We report data taken on 2 thin 
films grown at Rutgers University as reported elsewhere®. 
The data were taken approximately 30 minutes after the 
samples, originally sealed in vacuum, were received at 
JHU. The samples were mounted on a cold finger of a 
continuous He gas flow cryostat and cooled to 4 K in one 
hour. The thin films remained in He atmosphere for the 
remainder of the low temperature experiment. All data 
discussed in this letter were taken at a temperature below 
7 K. After the low temperature experiment, the samples 
were warmed up to room temperature in 5 hours, and 
then taken out of the cryostat and left exposed to atmo- 
spheric conditions for distinct time intervals, typically at 
least 24 hrs. 

The terahertz (THz) spectra were obtained in a home- 
built TDTS system at JHU between frequencies of 0.15 
to 2.50 THz. The THz beam was generated by excit- 
ing photo-carriers in a DC-biased Auston switch antenna 
with a pulsed infrared laser with pulses of approximately 
50 fs of width. The antenna emits radiation with a band- 
width of a few THz. The THz beam is collimated and 
then transmitted through the sample using a series of 
lenses and mirrors, and finally collected and focused onto 
another antenna where detection of the THz wave occurs. 
The detecting antenna is also excited by the pulsed in- 
frared beam and works analogously to the emitter, but 
in this case the photo-carriers are accelerated by the im- 
pinging THz wave rather than by a DC voltage. The 
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FIG. 1. (Color Online), (a) 
Real and (b) imaginary parts 
of the conductance of a thin 
film of 64 QL (f QL=9.4A) 
measured on different days af- 
ter growth, (c) and (d) same 
as (a) and (b) for a 128 QL 
film. It is clear from (a) and 
(c) that a dramatic decrease 
on the Drude spectral weight 
is observed as the film is fur- 
ther exposed to atmospheric 
conditions. 



current generated by the photo-carriers is proportional 
to the electric field of the THz beam. As the antenna 
only works while it is illuminated by the ultrashort in- 
frared pulse, one can sample the THz waveform in time 
by changing the relative time delay between the excit- 
ing pulses at the emitter and detector. This ability to 
measure the electric field as a function of time allows 
this technique to measure the full optical response, in 
the form of the complex conductivity, without the use of 
Kramcrs-Kronig transformation. 

Figure [l] shows the data from a 64 QL sample mea- 
sured over a span of 4 weeks, and for a 128 QL film 
measured for almost 3 weeks. Panel (a) shows the real 
part of the conductance for the 64 QL film, which is com- 
posed of a free-electron-like (Drude) response at the low- 
est frequencies and an absorption peak due to an opti- 
cal phonon around 2 THz. In our previous worrfES we 
found that, for samples that were kept at atmospheric 
conditions for at least 2 weeks, the Drude response did 
not depend on thickness, and we assigned it as a 2D 
response. For frequencies below the phonon peak, we 
find that if one tracks the conductivity of a newly grown 
sample with time, the conductance decreases with time 
as shown in panel (a) and (c). At the same time, it is 
clear that the scattering rate (half-width at half maxi- 
mum) of the Drude response is not affected significantly. 
This behavior is consistent with a decrease of the carrier 
density with time. The final value of the conductance 
at the last day of measurement, as well as its scattering 
rate, is essentially equal to the one reported before for the 



2D response in our study of the thickness dependence of 
the THz properties^!. We have shown previously that the 
spectral weight of this feature is completely accounted for 
by the expected carrier density of the topological surface 
states obtained from ARPESP^. Therefore, we assign the 
part that changes with time as the bulk response. Figure 
[ljb) shows the imaginary part of the conductance. This 
also shows clearly the Lorentzian line shape due to the 
phonon, as well as the free electron response as a pos- 
itive slope in the zero frequency limit. It is clear from 
these data that this slope decreases with time, indicating 
a loss of spectral weight with further exposure to ambi- 
ent conditions. Data for a 128 QL thin film shown in 
fig. [l|c) and (d) have the same qualitative behavior, and 
prove that the robustness of the 2D Drude response is an 
intrinsic feature of these topological insulator thin films. 

In order to quantify the changes in the spectra, we fit 
them with a model ac conductivity function that effec- 
tively parametrizes the different components in the re- 
sponse. There are 2 Drude terms, one representing the 
bulk (Dl) and the other the two surfaces of the film com- 
bined (D2), and there is a Drude-Lorentz (DL) oscillator 
for the response of the phonon: 

, v ^pDi/^ . ^ D2 /4tt i{ £oo -l)u 

(Txx{u) = ■ p + ■ p . 

ILU — 1 JJ1 tU! — 1 £)2 47T 

^DL - - iuT DL 
Here represents the high frequency optical transitions 
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FIG. 2. (Color Online), (a) and (b). Real part of the conductance of a thin film of 64 QL in days 1 and 25, respectively, as 
a function of frequency. The different contributions to the conductance are plotted as lines of different colors, (c) and (d). 
Imaginary part of the conductance of the same sample on the same days. 



contributions to the dielectric function. uj p Di and T d% are 
the Drude plasma frequency and scattering rate for the 
ith Drude component; u>dl, ^>pDL and Tjjl are, respec- 
tively, the center and plasma frequencies, and scattering 
rate of the optical phonon. Wc obtain the conductance 
by multiplying the ac conductivity by t the film's thick- 
ness, G xx — a xx x t. We show in figure [2ja) and (b) the 
results of the fits on the real part of the conductance with 
the different components plotted as individual traces for 
days I and 25, respectively. Fits for the other days are 
of similar quality, ft is clear that the largest change be- 
tween the 1st and 25th day of measurements occurs on 
the Drude component that we associate with the bulk 
response, and in particular it is the spectral weight (wj,), 
as opposed to the scattering rate, that is most affected. 
This Drude component corresponds to the bulk because 
we already know that the remaining contribution to the 
low frequency spectral weight after weeks of exposure be- 
haves two-dimensionally, as was shown previously^. In 
panels (c) and (d) we plot the imaginary conductances, 
where fits obtained simultaneously are of similar qual- 
ity. Therefore, the results of the fit show that the bulk 
response changes significantly with time as the film is 
exposed to atmospheric conditions, whereas the surface 
remains mostly unaffected. 

In Figure |3ja) and (b) we plot the time dependence 
of the fit parameters of the Drude components for the 
64 QL sample. After the first 2 days, the Drude spectral 
weight of the bulk in the 64 QL film decreases rapidly and 
after day 10 stays relatively constant. We note, however, 
that after day 10 the spectral weight in this Drude com- 
ponent is so small that it is very difficult to distinguish it 
from zero, therefore the error bars on the plot reach zero 
on these cases. The Drude component of the surface, on 
the other hand, decreases more gradually in the first few 
days, but then seems to level off for the remainder of 
the experiment, staying constant within the accuracy of 
the measurement. The scattering rate (proportional to 
the mobilities) of the surfaces seems to stay mostly un- 



changed within this time scale, indicating that the carrier 
density is the property that changes the most with ex- 
posure. We note that at the lowest measured frequencies 
and temperatures, for films that have aged the most, the 
surface contribution to the conductance can be up to 80% 
of the total conductance. This supports the scenario we 
proposed before^ that the THz response is dominated 
by the 2D carriers, and, surprisingly, allows us to con- 
clude that exposing the films to atmospheric conditions 
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FIG. 3. (Color Online), (a) Spectral weights (^) of bulk 
(black) and surface (red/gray) Drude components in a thin 
film of 64 QL as a function of measurement day. (b) Drude 
scattering rates for bulk and surface as a function of time for 
the same sample. 
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is a novel way of enhancing the surface to bulk ratio to 
the total film's conductance. 

The reduction of the bulk carrier concentration most 
probably results from the filling in of the Se vacan- 
cies present in the as-grown film by the oxygen in the 
atmosphere^. This effect would need to be stronger 
than the n-doping effect that water molecules have on 
the surface-dSI j n order to be consistent with our experi- 
ment. We would like to emphasize that the kind of aging 
effect we are reporting, which takes place over days in at- 
mosphere conditions, is clearly different than the kind of 
aging that takes place over hours in UHV*2U2I which has 
been shown to raise the surface Fermi level. The present 
case is probably driven by slow atomic diffusion into the 
bulk, while the UHV case is presumably driven by much 
faster physio- or chemisorption on the surface. We see 
no evidence for a third spectral feature that could derive 
from topologically trivial surface states. One might spec- 
ulate that these states, if they exist, have mobilities so 
low that they don't make an appreciable contribution to 



the conductance within our spectral window. 

In conclusion, we have shown that exposure of thin 
films of the topological insulator Bi 2 Se 3 to atmospheric 
conditions reduces the contribution of the bulk to the to- 
tal film's conductance by significantly changing the car- 
rier concentration. Using THz spectroscopy we were able 
to separate the effects on the mobility and carrier density 
of both the surface and bulk carriers, and thus conclude 
that in the time scales of days, for films exposed directly 
to the atmosphere, Se vacancies are replaced with Oxy- 
gen ions, which reduces the bulk carrier concentration 
without affecting the surface states properties that sig- 
nificantly. We think that the use THz spectroscopy will 
continue to aid in the study of the response of the topo- 
logical surface states and other 2D electron systems. 
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